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p ■ Abstract 

\ The Higgs boson is a central part of the electroweak theory and is crucial 

to generate masses for quarks, leptons and the weak gauge bosons. We use 
a 4-dimensional Euclidean lattice formulation of the Higgs- Yukawa sector of 
the electroweak model to compute physical quantities in the path integral 
approach which is evaluated by means of Monte Carlo simulations thus al- 
lowing for fully non perturbative calculations. The chiral symmetry of the 
model is incorporated by using the Neuberger overlap Dirac operator. The 
here considered Higgs- Yukawa model does not involve the weak gauge bosons 
and furthermore, only a degenerate doublet of top- and bottom quarks are 
incorporated. The goal of this work is to study the resonance properties of 
the Higgs boson and its sensitivity to the strength of the quartic self coupling. 
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1. Introduction 

The model under consideration is the pure Higgs- Yukawa sector of the 
electroweak standard model formulated on a 4-dimensional Euclidean space- 
time lattice. The field content of the model consists of a complex scalar 
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Higgs doublet and a fermion doublet coupled in a chirally invariant way to 
the scalar sector. Our lattice model closely resembles the weak interaction 
of the continuum standard model in the limit where the gauge fields are 
switched off. The chiral nature of the weak interaction is realized in our 
lattice setup with the help of the Neuberger overlap Dirac operator [l| which 
allows to incorporate an exact chiral symmetry on the lattice 

Within the past years renewed efforts have been undertaken to investigate 
the lower and upper Higgs boson mass bounds within this model in the 
framework of lattice field theory 0, 0, Isj , see these references also for accounts 
of earlier lattice investigations of Higgs boson mass bounds. However, the 
decay width of the Higgs boson was not taken into account in refs. [1, 0] 
assuming that its effect on the Higgs boson mass is small. The first main 
aim of this work is to treat the Higgs boson as a true resonance and to 
investigate the unstable nature of the Higgs boson from first principles. In 
this way, the assumption of the preceding work js], 0] can be tested and it 
can be investigated, whether the Higgs boson mass bounds are affected. The 
second goal is to determine the Higgs boson width as a function of the quartic 
self-interaction of the Higgs fields, including large values of the coupling. In 
particular, we want to investigate, whether in the decay channel of the Higgs 
boson into massive vector bosons the Higgs boson becomes a broad resonance 
at large couplings, or whether the width is still rather narrow. 

The strategy to determine the Higgs boson resonance parameters is the 
analysis of the scattering phases. The computation of the scattering phase 
within lattice field theory is in turn based on the determination of the vol- 
ume dependence of energy eigenvalues |6|. An extension of the finite size 
method of ref. [oj, which was formulated in the center of mass frame, has 
been proposed in [7| and is based on the analysis of the energy levels within 
a moving frame. The method is complementary to the original method and 
allows to obtain more energy levels on the same underlying field configura- 
tions. We employ this method here to reduce the errors on our extraction 
of the scattering width and the resonance mass, which we are then able to 
determine at the 0(10%) level. 

2. The model 

The Higgs- Yukawa model is defined by the Lagrangian and the corre- 
sponding generating functional for the Green functions of the theory. With 
regard to the lattice formulation of the model, the Euclidean version of the 
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model will be considered here. The particle content contains the scalar sec- 
tor and the heaviest quark doublet consisting of the top and the bottom 
quark. Due to the neglection of gauge bosons the Lagrangian exhibits a 
global SUwC^) x ^y(l) inner symmetry rather than a local (gauge) symme- 
try. The Euclidean action is given by 



+ tlpt + blpb + yt,(^^-ipbn + yt(^^-^tR +h.c.. (1) 



L \ / L 

i2 



Here, m , A and yt(b) denote the bare scaler mass, the bare quartic coupling 
and the bare top (bottom) Yukawa coupling, respectively. It is common in 
lattice field theory to rewrite the scalar sector by rescaling the fields with 
a factor a/2k and furthermore, the scalar doublet can be expressed as a 
quarternion 



'^2 ^2 



\-%a^4P, jG{1...3}, 0''gM. (2) 



The rescaled fields are: 

$^ := ^= d>^. 

The scalar fields in the usual notation of eq. ([T]) can be recovered by identi- 
fying 



where V2i,V'2 are the original, coniplex valued scalar fields. 

In 1982 Ginsparg and Wilson j8| proposed a relation which defines a class 
of lattice Dirac operators and which is since then known as the Ginsparg- 
Wilson relation 

75!) + 1^75 = aD^^RD. (3) 

Here a is the lattice spacing and i? is a positive constant. The Ginsparg- 
Wilson relation can be utilized to construct a lattice modified chiral symme- 
try which recovers the desired continuum chiral symmetry in the limit a — )■ 
[2|. The modified lattice projectors for the left and right handed components 
are given by 

P±:=^(l±75), 75 := 75(1- ai?D). 
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The modified projectors are used to define the chiral components of the 
spinor fields on the lattice 



%Ijr = P+ij, IpL = P-lp, 

As in the continuum theory the free part of the fermion Lagrangian can 
be written as a sum of the left and right handed lattice spinor fields 

The modified projectors allow to formulate a theory of continuum-like left- 
and right handed fermions on the lattice for any non-vanishing value of the 
lattice spacing a. 

In order to specify the lattice action, a Ginsparg- Wilson type Dirac op- 
erator has to be introduced. The here presented results are based on the 
Neuberger overlap operator which satisfies the Ginsparg- Wilson relation 
and is given by 



pM = ^/l + ^Ll, A = D^-^, 0<p<2r 
a y y/A'fA) a 



(4) 



p is chosen to be 1 in this work and denotes the Wilson Dirac operator. In 
it has been demonstrated that the Neuberger overlap operator is, despite 

its complicated form, a local operator. 

The full Euclidean discretized action is given by 



y 

x,y ) \by 



(5) 



In the following, the lattice spacing a is set to unity. It can be given 
a physical value by identifying the vacuum expection value on the lattice, 
afphys with the measured value of fphys = 246 GeV. As mentioned before, 
the couplings are scaled by the k parameter. The parameterization of the 
Lagrangian in ([T]) can by recovered with the identities: 

A 2 l-2NfX-8K yt,i> 

^ = ^0 = ' yt,b - 
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The here considered Higgs- Yukawa model has aheady been studied ex- 



tensively. In refs. [10|, lUl the phase structure has been determined, lower 
and upper Higgs boson mass bounds have been obtained in 0, 0] and also an 
extension to a possible fourth quark generation has been discussed in |l2 . 
Many details of the model, _its numerical implementation and algortihmic 
aspects can be found in ref. 
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3. Resonances from lattice simulations 



The method to compute the resonance parameters in finite volume was 
proposed in [6|. It has been successfully employed for the case of the Higgs 
boson mass in the pure theory IJ]; here we want to extend the method 
to the Higgs- Yukawa model as considered in this work. 

Let us start with the optical theorem which states that the total cross 
section in the elastic region of forward scattering is given by 



(6) 



j=0 



q denotes the centre of mass momentum and s is the centre of mass energy of 
the two incomming particles with four momentum pa and pb, s = {pa+ Phf ■ 
T/i is the spherical decomposition of the interacting part T of the S-Matrix 

S = I + iT. 

The total cross section exhibits a peak near a resonance which can be well 
parameterized with a Breit-Wigner function containing the resonance mass 
m and the width F, 



p2 _ jyip. _|_ j^^Y 



i=o 



l)2sin^5j(s) 



(7) 



where we have also expressed the total cross section in eq. (E]) in terms of the 
scattering phases 5j. 

In order to develop a general method to extract the scattering phase from 
the energy levels in a finite box, it is sufficient to investigate the problem in 
non-relativistic quantum mechanics. The non-relativistic result can then be 
transferred to the case of quantum field theory. This remarkable result has 



been demonstrated in refs. 15, 16 



5 



As a result of these works the computation of the scattering phase Sq in 
which we are interested in, proceeds as follows. The first step is to compute 
the energies of the two-particle massive Goldstone system, W = 2a/ + /c^, 
from which the value of k can be extracted. Knowing k, the scattering phase 
(5o is obtained from 

Zoo (1; q^) 27r 

Zoo = Y] -7i=— 2- 

The arialysis of scattering phases has been extended to moving frames 
in 0, [13, 18|, where one of the two particles is at rest. The method is 



complementary to the centre of mass frame and allows to compute more 
data for the scattering phases from the same configurations. The obtained 
two particle energies have to be translated back to the centre of mass frame. 
In general, the choice of a moving frame implies that one has to consider 
irreducible representations of a subgroup of the cubic group. The remaining 
symmetry depends on the selected directions of the moving frame. The 
here used irreducible representations are given in [Tj]. The modification of 
the relation between the two particle energy in the moving frame and the 
scattering phase is given by 

tan^o g = ^JTT^' ^ = ^ ' 

where p* denotes the momentum which has been transferred back to the 
centre of mass frame by a Lorentz boost. The modified Zeta function is 
defined by 



2 _ /,2^2 ' 



4vr (r2 - q 
Prf := <J f G M^|7"^ (n + ^d) ,neZ^ 

The vector d is related to the total momentum of the moving frame P. 7 
is the usual Lorentz factor. Below are some definitions which are needed to 
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compute the modified Zeta function 



P 





P 
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Wl denotes tlie two particle energy in tlie moving frame wliich can be com- 
puted with time shce correlators. The corresponding observables are defined 
in the next section. n\\ and n_L is a decomposition of the vector n in its 
parallel and perpendicular parts with respect to the centre of mass velocity 

V 



4. Numerical results 

The Higgs boson decays dominantly to any even number of Goldstone 
bosons, if kinematically allowed. The physical set-up chosen here allows al- 
ways for such a decay, i.e. we were always working in the elastic scattering 
region. In order to achieve this situation, external currents were introduced, 
providing appropriate masses to the else mass-less Goldstone bosons. The 
bare parameters used in our simulations are given in table [TJ The corre- 
sponding physical value of the cut-off, the Higgs boson propagator mass, the 
Goldstone boson mass and the obtained top quark mass are summarized in 
table El 

The Goldstone theorem ensures that the Goldstone bosons are massless. 
Due to an external current J which couples to one of the components of 
the scalar fields in the complex SUw{'^) doublet, the symmetry is broken 
explicitly in the Lagrangian. The Goldstone bosons acquire a mass and they 
form a vector under cubic rotations. The magnitude of the current J is 
chosen such that the ratio of the Higgs boson mass to the Goldstone 
boson mass Mq is roughly 3. Here and below the superscript p in and 
Mq denotes that the mass was extracted from the analysis of the momentum 
space propagator. This numerically computed propagator was then fitted 
to a formula motivated from perturbation theory and the real part of the 



n ■ v)v 



n± := n — n\\. 
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Table 1: The table summarizes the bare parameters for the Monte Carlo sim- 
ulations which were performed in order to determine the scattering phases. 
J denotes an external current which is coupled to one component of the 
scalar field. J is adjusted such that the the ratio of the Higgs boson mass 
and and the Goldstone boson mass Mq is approximately 3 (see ta- 
ble [2|), which ensures that the Higgs boson decays into any even number of 
Goldstone bosons. 



A 




y 


J 


0.01 


0.1295 


0.36274 


0.001 


1.0 


0.2445 


0.49798 


0.002 


oo 


0.3020 


0.57390 


0.002 



complex pole of this fit function has been identified with the Higgs boson 
mass, see refs. js], 0, 12, 13 . 



The Higgs field is a singlet under cubic rotations and transforms as ele- 
ments in the representation. The two particle energies discussed in the 
previous section are constructed from the two particle Goldstone singlet as 
it has the same quantum numbers as the Higgs boson. 

The analysis of the resonance parameters involves several lattice volumes 
with identical bare parameters in order to compute the momentum depen- 
dence of the scattering phase. As shown in table [H there are three distinct 
set of simulation parameters which shall be characterized with the value of 
the bare quartic coupling (A G {0.01, 1.0, oo}). For each of the three values 
of the quartic coupling the simulations were performed on lattice volumes up 
to 40^. In the following the two particle energies of the two Goldstone boson 
states will be discussed. Once these energy levels are known, the unstable 
nature of the Higgs boson can be studied by the method described in the 
preceeding chapter. 

The Goldstone bosons are stable particles such that their ground state 
energy can be calculated from the two point time correlation function. The 
concept of time correlators is widely used in lattice field theory and there 
are reliable techniques to extract mass eigenvalues from such correlators. 
The method of choice in this work, is the analysis of the correlation matrix 



15L Il9| . The correlation matrix built from an operator Oaif) is given by 

Ca^m ■■=T (o.(t)o,(t'))e , 



t-t' =At 
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Table 2: The table summarizes the obtained masses in lattice units for the 
bare parameters given in table [TJ The Higgs and Goldstone boson masses, 
and respectively, are extracted from the propagator neglecting 
the width. The quark mass is computed from the corresponding time slice 
correlator and the last columns show the renormalized vacuum expectation 
value vr and the cut-off (A). For the smaller lattices we have typically 
used 20000 and for the largest lattice typically 10000 configurations. The 
autocorrelation times were about one in all cases, see also [1.3]. Note that, 
due to a technical problem, for the determination of Mfj at A = 1.0 we 
have used a significantly smaller statistics than at the other simulation 
points. 



A 






mt [GeV] 


avR 


A [GeV] 


0.01 


0.278(1) 


0.085(2) 


174(1) 


0.2786(3) 


883(1) 


1.0 


0.386(28) 


0.133(4) 


179(2) 


0.1637(5) 


1503(5) 


oo 


0.405(4) 


0.129(1) 


178(1) 


0.1539(2) 


1598(2) 



where Lt is the temporal size of the lattice. Throughout this chapter the 
temporal extent will be Lt = 40. The subscript c denotes that the discon- 



nected part of the correlator has been subtracted. It has been shown in [15 
that the eigenvalues of the correlation matrix decay exponentially with rising 
time separation At. In the following the operators which contribute to the 
two Goldstone system are collected. 

The definition of the observables in the centre of mass frame is given by 

Oo{t) := H{0,t) 

Oi(t):=i=^5^^^(n,t)^(-n,t) 

fieQi 

Qi = {n G = 0} , |Qi| = l 

Q2 = {neZ^\n^ = 1} , \Q2\=Q 
1 1 



Q^ = {neZ''\n^ = 2] , IQsl = 12. 



9 



The correlation matrix is thus a 4 x 4 matrix. H{0, t) and Q{n, t) denote the 
standard, zero momentum projected Higgs and Goldstone boson interpolat- 



ing fields H{x), Q{x), see e.g. |l4j. Il3 



In order to collect more data for the scattering phases, the modification 
of the method to a moving frame was analysed as well. The moving frame is 
characterised by a constant vector d which indicates the momentum of the 
frame. The observables for the moving frame are constructed such that one 
of the Goldstone bosons is at rest while the other can take any momentum 
allowed on the lattice. The selection of a constant vector d breaks the cubic 
symmetry and thus special care is needed while constructing the observables. 
Fortunately, it turns out that the Af sector does not need much modifica- 
tion and explicit relations are given in [3]. The lowest energy eigenstates 
are associated to the lowest possible relative momentum and thus only mov- 
ing frames with momentum d = (0, 0, 1) and permutations thereof will be 
considered. The observables are (cj is the unit three-vector in direction i) 

di €-i 

0,^Jt) ■.= H{lt) 
0,^Jt) :=^^(ti:,t)^(0,t) 

0,-;_2(At) := i Yl Q^in + d,,t)g{-n,t) 



Qd„2 = 1^ G Z|n ■ (i* = 0, ra^ = l| 



^d-.sW ■=g''{2di,t)g{-di,t). 

The correlation matrix in the moving frame is then given by 

* \t-t'\=At i=l 

a,/3e {1,2,3}. 

The energy levels obtained from the moving frame are connected to energy 
levels in the corresponding centre of mass frame by Lorentz transformation. 

Fig. [Ushows the obtained cross sections and Fig. |2]the obtained scattering 
phases for the three different physical situations (different values of A) we 
have used. The cross sections are plotted against the energy Wk while the 
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scattering phase, which take values in the interval [0,7r], are plotted against 
the momentum k. If the scattering phase 6{k) passes through | it indicates 
the existence of a resonance. Hence, all three setups involve an unstable 
Higgs boson. The cross section can be decomposed into spherical harmonics 
and is parametrised by the scattering phase 

4-7r 

^(^)= pE(2j + l)sin^(^.(^)) (8) 

j=0 

which becomes 

a{k) ~ -^sin^ ('^o(fc)), 

when the contribution of the higher angular momenta j > are neglected. 
We use a fit function to describe our data for a{k) which is motivated by a 
Breit-Wigner form and allows us to determine the resonance mass Mh and 
the width Vh'- 

f{k) := IGtt j^^^ ^. (9) 

{Ml - Ami) [m - Mlf + MlTl) 

In fig. [H the solid curve shows the result of the fit for the cross section 
and in fig. |2] we show the corresponding scattering phases. The scattering 
phase is directly related to the cross section in equation ([9]) as can be seen in 
equation ([7]). Given the scattering phases, the cross section can be obtained 
from equation ([H]) setting j = 0. 

Finally table E] summarizes the results for the Higgs boson mass obtained 
by different approaches. The physical (resonance) Higgs boson mass aMn is 
compared to the mass obtained from the Higgs boson propagator aM^ and 
the energy eigenvalue from the correlation matrix analysis (GEVP) which 
corresponds to the Higgs boson mass. The results for aM^ and the GEVP 
were obtained after an extrapolation to infinite volume, see fig. [3] for the ex- 
ample of the infinite volume extrapolation of the eigenvalue from the GEVP. 
For completeness, we also show the finite size behaviour of the renormalized 
vacuum expectation value and the top quark mass in fig. HI 

The simulations at A = 1.0 and A = oo belong to a cut-off of around 1.5 
TeV. It was necessary to reduce the cut off to 880 GeV for the smallest quartic 
coupling A = 0.01 in order to meet the resonance condition {-^ ^3). As 
one can see from table [31 the Higgs boson mass is still well below the cut off. 
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0.32 0.36 


0.4 0.44 0.4S 0.52 

(b) 




0.32 0.3ti (1,4 0,44 0.48 



(c) 



Figure 1: The figure shows the cross-sections as obtained in three different physical 
situations for A = 0.01 (a), A = 1.0 (b) and A = oo (c). The cross sections are 
plotted against the energy Wfc. The red points refer to cross sections obtained 
from the analysis in the centre of mass frame as originally proposed in |6|. The 
blue points denote the cross sections computed within a moving frame [7j. The 
solid line is a fit to the data using eq. ([9]). The computations were performed on 
various lattice volumes x 40 where G {12, 16, 18, 20, 24, 32, 40} . 



Especially for the smallest quartic coupling A = 0.01 the resonance region is 
very small (mc = 0.09(1) ^ 0.18 < Wk < 0.36) which in turn necessitates 
large lattice volumes in order to obtain energy eigenvalues which lead to 
scattering phases near the resonance mass. The plots in fig. [1] and fig. [2] 
show that the analysis of the moving frame can help significantly to extract 
reliable results. 
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0.06 0.08 0.1 0.12 0.14 0.16 0.04 0.08 0.12 fl.lfi 0.2 0.24 0.04 0.08 0.12 O.lfi 0.2 0.24 

k k k 



(a) (b) (c) 

Figure 2: The figure shows the scattering phases that correspond to the cross 
sections shown in fig.[T]from which we also take the notation. The scattering phases 
are plotted as a function of the momentum k. The vertical dotted line indicates 
the inelastic threshold. The solid line shows the scattering phase obtained from 
the fit function to the cross sections displayed in fig. [TJ 



Table 3: The table summarizes the obtained final results on the resonance 
mass Mh and the resonance width Th of the Higgs boson. A denotes the 
bare quartic coupling. A is the cut-off of the theory. The following two 
columns display the resonance parameters computed from the scattering 
phases. The P^r degree of freedom from the fits to obtain the resonance 
parameters are x^/d.o.f. = 1.3, x^/d.o.f. = 1.0 and ^jd.o.f. = 1.2 for 
A = 0.01, A = 1.0 and A = oo, respectively. F^'^* is the width obtained from 
perturbation theory where a non vanishing mass for the Goldstone bosons 
has been considered (2o|. Finally the mass extracted from the propagator 
as well as the eigenvalue corresponding to the Higgs boson mass computed 
with the help of the correlation matrix is shown. 



A 


A [GEV] 










GEVP 


0.01 
1.0 

oo 


883(1) 
1503(5) 
1598(2) 


0.278(3) 
0.383(6) 
0.403(6) 


0.0018(14) 
0.0169(4) 
0.037(9) 


0.0054(1) 
0.036(8) 
0.052(2) 


0.278(2) 
0.386(28) 
0.405(4) 


0.274(4) 
0.372(4) 
0.403(7) 
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0.32 - 




0.26 I 1 1 1 1 

0.002 0.004 0.006 0.008 



Figure 3: The figure shows the energy eigenvalues corresponding to the Higgs bo- 
son mass obtained by the correlation matrix analysis. The eigenvalues are plotted 
against the inverse squared spatial lattice extent in order to perform an extrapola- 
tion to infinite volume. From top to bottom the data correspond to A = oo, A = 1 
and A = 0.01 



0.3 
0.28 
0.26 I- 
0.24 . 
0.22 

0.2 
0.18 - 
0.16 
0.14 
0.12 





0.004 
l/i^ 



0.16 
0.14 
0.12 
0.1 
0.08 



(a) 



1/L= 

(b) 



Figure 4: The figure shows the finite size effects of the scalar vacuum expectation 
value and the fermion mass. The top quark mass is computed from the fermion 
time slice correlation function. The figure shows also an extrapolation to infinite 
volume starting from lattice volumes of at least 16^. From top to bottom the data 
correspond to A = 0.01, A = 1 and A = oo 
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5. Conclusion 



In this letter we have investigated the resonance parameters of the stan- 
dard model Higgs boson from simulations on a Euclidean lattice in the limit 
of the electroweak theory where the gauge fields are neglected. In this situa- 
tion the Higgs boson decays into massive vector bosons. The resonance mass 
and the width of the Higgs boson is computed using the finite size method 
proposed in [6i], applicable for the center of mass frame, combined with the 
extension to a moving frame as proposed in Using both kind of Lorentz 
frames a large number of scattering phases could be obtained which allowed 
us to compute the Higgs boson resonance mass and width with an accuracy 
of 0(1%) for the resonance mass and 0(10%) for the width, see table Ei 

We have employed three values of the quartic coupling A of the Higgs 
field ranging from small perturbative values to A = oo with parameters cor- 
responding to a large renormalized quartic coupling. In all three cases the 
Higgs boson width is not larger than about 10% with respect to the reso- 
nance mass. Therefore, the corresponding total cross section exhibits a clear 
resonance peak even at the strongest value of the quartic coupling. 

The values of the resonance mass we have extracted here are in very good 
agreement with earlier determinations js], 0] where the Higgs boson width has 
been neglected. This finding provides confidence to and justifies a posteriori 
the Higgs boson mass bounds determined in js], 0]- In addition, a comparison 
of the Higgs boson width with results from perturbation theory reveals a very 
good agreement as can be seen in table |3l 

Figure [5] summarizes the obtained total cross sections for the three dif- 
ferent values of the bare quartic couplings. The displayed curves correspond 
to the fits shown in figure [T] using the parameterization of eq. It will be 
interesting to apply the techniques used in this paper for the investigation of 
the Higgs boson mass width in presence of a possible fourth fermion gener- 
ation, where non-perturbative effects might appear for very heavy masses of 
the four th g eneration quarks. Simulations in this direction have been started 



already [12 
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k 

Figure 5: The figure shows the total cross section of the Higgs boson within the Higgs- 
Yukawa model as a function of the momentum k. With regard to the standard model, 
this cross section is associated to the decay of the Higgs boson into the weak gauge bosons 
W^,Z. The highest peak belongs to the smallest bare quartic coupling A — 0.01 and 
corresponds to a Higgs boson resonance mass of Mh — 248 ± 1 GeV and a resonance 
width of Tjf = 6.2 ± 0.9 GeV. The next peak is obtained at A = 1.0 and corresponds to 
Mh — 562 ± 2 GeV and Th = 50 ± 6 GeV. The last peak is associated to infinite bare 
quartic coupling and corresponds to Mh = 618 ± 5 GeV and Th = 60 ± 6 GeV. In all 
cases, the resonance width is less than 10% with respect to the resonance mass and thus 
the corresponding total cross section exhibits a clear resonance structure. The energy 
ranges corresponding to the different A values can be read off from figure [TJ 
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